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ABSTRACT 


The  finite  difference  (FD)  method  yields  solutions  to  discretized  versions  of  the  full  acoustic  wave  equation  for 
arbitrarily  complex  media.  The  method  is  reliable  at  all  angles  of  propagation,  including  backscatter.  This  offers  an 
advantage  over  other  standard  propagation  methods  in  wide  use,  as  it  allows  for  accurate  computation  of  acoustic 
energy  levels  in  the  case  where  significant  scattering  can  occur  near  the  source,  such  as  may  happen  for  an  explosion 
near  the  surface,  or  underground.  It  also  allows  for  the  investigation  of  the  penetration  of  infrasound  energy  into 
classical  shadow  zones,  where  ray  theory  predicts  the  upward  refraction  of  sound.  This  fits  in  with  nuclear 
monitoring  goals,  in  that  it  allows  for  an  improved  understanding  of  the  generation  and  propagation  of  infrasound 
energy  from  arbitrary  sources,  including  underground  and  near- surface  explosions. 

The  effects  of  1)  wind,  2)  attenuation,  and  3)  gravity  on  infrasound  propagation  are  examined  separately,  using 
finite-difference  time-domain  methods.  The  method  is  applied  in  2-D  Cartesian  coordinates.  The  method,  including 
the  effects  of  wind,  is  applied  to  realistic  problems  in  infrasound  propagation,  including  the  penetration  of 
infrasound  energy  into  “shadow”  zones,  as  defined  by  ray  theory.  The  effects  of  diffraction  and  topography  are 
examined  for  sources  near  the  surface.  It  is  shown  that  the  FD  approach  can  be  used  to  solve  for  sound  intensities  in 
complex  models  that  may  include  high  material  contrasts  and  arbitrary  topography. 
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OBJECTIVES 


A  finite  difference  time  domain  (FDTD)  method  is  sought  for  the  numerically  propagation  of  infrasound  energy 
through  realistic  atmosphere  properties,  including  variable  winds,  sound  speeds,  and  density  variations,  and  the 
interaction  of  the  sound  field  with  topography.  The  FDTD  method  yields  the  solution  to  a  discretized  version  of  the 
full  acoustic  wave  equation  for  arbitrarily  complex  media,  and  can  accurately  handle  propagation  reliable  at  all 
angles  of  propagation,  including  backscatter.  This  offers  an  advantage  over  other  standard  propagation  methods  in 
wide  use,  as  it  allows  for  an  accurate  computation  in  cases  where  significant  scattering  can  occur  near  the  source, 
such  as  may  happen  for  an  explosion  near  the  surface  or  underground.  The  effects  of  wind  or  gravity  on  infrasound 
propagation  may  also  be  incorporated  into  the  infrasound  propagation  problem  with  relative  ease  using  FD 
techniques.  In  this  paper,  a  method  to  handle  the  inclusion  of  attenuation  into  the  time-domain  computations  is 
investigated,  as  well  as  the  effects  of  topography  on  infrasound. 

The  equation  governing  the  effects  of  gravity,  wind,  and  attenuation  on  the  infrasound  field  are  derived  in  this  paper, 
and  their  effects  on  propagation  through  the  atmosphere  are  demonstrated.  The  penetration  of  infrasound  into  a 
classical  “shadow  zone”  is  also  demonstrated  both  for  a  model  of  a  source  above  a  flat  surface  with  high  impedance 
contrast,  and  also  for  the  same  source,  above  a  rough  surface  with  high  impedance  contrast.  In  this  way,  we  examine 
the  effects  of  topography  on  propagation  of  infrasound  energy  into  a  shadow  zone. 

RESEARCH  ACCOMPLISHED 


Equations  governing  infrasound  propagation  through  the  atmosphere  are  derived  here.  The  effects  of  gravity,  wind, 
and  atmospheric  attenuation  of  the  acoustic  field  are  examined  separately.  Following  that,  the  effects  of  topography 
on  infrasound  penetration  into  an  acoustic  shadow  zone  are  shown.  In  each  case,  the  methods  are  applied  in  2-D 
Cartesian  coordinates. 


Linear  acoustic  propagation  equations 


The  equations  governing  propagation  of  sound  in  the  atmosphere  are 

i)  the  conservation  of  momentum, 


Dv 

P—  =  -Vp  +  F,  (1) 

Dt 

which  relates  the  particle  velocity,  V ,  to  the  total  pressure,  p,  and  density,  p,  and  the  external  volume  force  per  unit 
mass,  F ; 

ii)  the  conservation  of  mass, 


iii)  and  the  equation  of  state, 


Dp 

Dt 


-pV  •  v, 


(2) 


Dp  „2  DP 


(3) 


Dt  Dt 

where  c  is  the  total  sound  speed  in  the  presence  of  acoustic  perturbations.  It’s  assumed  in  the  equation  above  that 
infrasound  propagation  is  isentropic,  i.e.,  that  entropy  is  constant.  That  is,  it’s  assumed  that  attenuation  is 
sufficiently  small  so  that  wave  propagation  dominates  over  diffusive  propagation.  This  is  a  valid  approximation  for 
infrasound  propagation  within  a  realistic  atmosphere. 


The  convective  derivative  (also  known  as  the  Lagrangian  derivative)  acting  on  a  quantity  G  is  defined  as 

DG  dG 

- = - +  (v«  V)G, 

Dt  dt 


(4) 


where  the  derivative  on  the  left  represents  the  change  of  G  with  time  in  a  frame  moving  with  the  fluid.  The  first  term 
on  the  right  side  of  Equation  4  represents  the  change  in  G  at  a  point  fixed  in  space.  The  second  term,  called  the 
advective  term,  represents  the  change  in  G  as  the  observer  moves  with  the  fluid  at  the  velocity  v  .  For  infrasound 
investigations,  quantities  are  expressed  in  terms  of  a  fixed  point  in  space  in  order  to  compare  computational  results 
with  observations  made  at  stationary  sensors. 
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Propagation  of  sound  waves  in  the  atmosphere  introduces  fluctuations  in  the  pressure,  density,  and  particle  velocity 
values.  The  standard  procedure  in  solving  Eqs.  1  through  3  is  to  consider  a  perturbed  solution  of  the  form 

2  2  2 

P  =  Po+Ps>  p  =  Po  +  Ps’  V  =  V0  +  VS’  c  =co+cs ;  (5) 

where  pQ,  p0,  and  vQ  are  the  ambient  values  of  pressure,  density,  and  particle  velocity  in  the  absence  of 
perturbations,  and  ps,  ps ,  and  are  perturbations  caused  by  the  passage  of  a  sound  wave.  The  pressure  fluctuations 

cause  slight  fluctuations  in  the  squared  sound  speed.  Waveforms  are  derived  by  computing  the  pressure 
perturbations,  ps,  as  a  function  of  time.  Linear  expressions  for  the  particle  velocity,  vs,  and  pressure  and  density 
fluctuations,  ps  and  ps ,  are  derived  by  retaining  only  first  order  terms  in  the  pressure,  velocity,  and  density.  The 
zeroth  order  terms  cancel  because  the  ambient  atmosphere  is  assumed  to  be  in  equilibrium. 


Linearized  acoustic  propagation ,  including  gravitational  effects 


In  this  sub-section,  we  consider  solutions  in  a  static  medium,  thus  V0,  denoting  the  wind  velocity,  is  equal  to  zero.  If 
the  effects  of  gravity  are  included  in  the  propagation  equations,  the  external  volume  force  per  unit  mass  may  be 
expressed  as 

F  ~  ~(Po  +  Ps)g  +  /’  (6) 

where  the  gravity  g=  [0,0,9. 8m/s2],  and  f  represents  all  external  forces  other  than  gravity.  Using  Eqs.  5  and  6  in 
Eq.  1 ,  and  linearizing,  yields 

dv* 

P  0  —  =  -Vps-Psg  +  f-  (7) 

dt 

under  the  assumption  that  there  is  no  wind,  and  using  the  hydrostatic  equation  VpQ  =  -pQg .  Combining  Eqs.  2  and 
3  yields  the  following  expression  for  pressure  perturbations 


fys 

dt 


=  Povs*S-Pocov> 


(8) 


where  the  hydrostatic  equation  has  again  been  used.  Horizontal  gradients  in  the  ambient  pressure  and  density,  pQ  and 
p0 ,  are  negligible  in  the  horizontal  direction.  Equation  2  becomes 


5P.y 

dt 


-  + vc 


vPc 


l 


dp 


L  dt 


S  +  vs  •  yPo 


(9) 


Combining  the  ideal  gas  law  (pQ  =  PoRT)  with  the  hydrostatic  equation,  and  the  relation  between  sound  speed  and 
2 

temperature  (  cQ  =y. RT )  yields 


dpg_ 

dz 


YP 0S  2p o  dc 


cQ  dz 


(10) 


where  y  «  1.402,  T  is  the  temperature  in  degrees  Kelvin  and  R  is  the  gas  constant  for  dry  air,  R  =  287.04  J  kg'1  K'1. 
Using  the  hydrostatic  equation  for  vertical  pressure  gradients  and  Eq.  10  for  vertical  density  gradients,  Eq.  9  can  be 
expressed  as 


dps  1  8ps 


dt 


C20  dt 


+ 


PoVS' 


(y  -  !)£  ,  2 

- —  +  —  Vco 

-  ct 


(11) 


~o  " o 

Equations  7,  8,  and  11  are  a  complete  set  of  equations,  in  a  form  amenable  to  FD  time  domain  computations. 


Note  that  Eq.  1 1  may  also  be  expressed  as  (Gill,  1982) 


Sp5  ^  _  ,t2  1  ops 


dt 


-PowsN  /S  = 


cl  St 


(12) 


where  the  Brunt- Vaisala  frequency  N — or  buoyancy  frequency — is  defined  by 
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N2=g 


(r  -  i)g. 


2 

+  —  Vc 
c  O 
o 


(13) 


J 


The  buoyancy  frequency  N  is  on  the  order  of  1/300  Hz  within  the  atmosphere.  The  equations  suggest  that  the  effect 
of  buoyancy  is  negligible  at  the  upper  end  of  the  infrasound  range  but  becomes  more  significant  with  decreasing 
frequency.  Results  of  computations  are  compared  below  for  sources  with  peak  frequencies  of  0.2  Hz  and  0.015  Hz, 
for  gz=  0  and  gz  =9.8m2/s2.  In  each  case  the  sound  speeds  are  assumed  to  be  uniform  (i.e.,  an  isothermal  atmosphere) 
and  exponentially  decreasing  density.  The  results  confirm  that  gravitational  effects  are  insignificant  at  higher 
frequency. 


Figure  1.  Results  shown  above  are  for  models  with  exponentially  decreasing  density,  a)  The  source  peak 
frequency  is  0.2  Hz.  There  is  less  than  0.2%  difference  between  the  solution  for  g=0  and  for 
g=9.8m2/s2.  B)  The  case  for  a  source  peak  frequency  of  0.015  Hz.  There  is  less  than  3.5%  difference 
between  the  solution  for  g=0  and  for  g=9.8m2/s2. 


Linearized  acoustic  propagation ,  including  wind  and  gravitational  effects 


When  wind  is  included,  i.e.,  v0  ^  0,  we  use  Eq.  4,  the  convective  derivative  to  find 

dv«  ____  _ 

Po(—  +  vs  •  Vvo  +  vo  •  Vvs)  +  P.sv0  •  Vvo  =  ~vPs  -  P^g  +  />  (14) 

ot 

for  particle  velocity.  Re-arranging,  and  assuming  that  wind  is  invariant  along  path  (  v0  •  Vv0  =  0 ,  which  is  true  for 
negligible  turbulence)  this  becomes 

.  1  - 

—  =  -vs  •  Vv0  -  v0  •  Vvs - (S7ps  +  P5g  -  /),  (15) 

8t  pQ 

i.e.  wind  shear  is  included  in  this  formulation.  The  above  result  differs  from  Ostashev  et.al.  (2005)  in  that  gravity  is 
included  and  thus  density  perturbations  must  be  computed. 


The  equation  for  pressure  perturbations  becomes 

dp, 


—  +  vs  *  vPo  +  vo  *  vPs 

Ot 


2  2  2 
-P oco  v  •  vs  -  P ocs  v  •  vo  -  P sco  v  •  vc 


(16) 


It’s  assumed  that  turbulence  is  insignificant  (as  in  Ostashev  et.al.,  2005)  so  V  •  v0  is  negligible.  The  hydrostatic 
equation  is  used  to  get  the  following  form  for  pressure  perturbations: 


dPs  z 

—  =  PoV5,g-V0*  V^-P0C0V#VS 
Ot 


(17) 


The  equation  for  density  perturbations  becomes 


dp  s 
.  dt 


+  VS*  VPo+V0,VP5 


2 

+  c5(vo*  vPo)  =  —  +  vs' 
ot 


VPo  +  v0  •  Vps 


(18) 
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Winds  are  horizontal  and  the  ambient  density  varies  only  with  altitude,  so  v0  •  Vp0  =  0 .  So,  using  the  hydrostatic 
equation  for  Vp0  and  Equation  10  for  the  spatial  density  derivative  Vp^and  rearranging  yields 


dt 


1 


dp„ 

—  +  v0  •  Vps 
_  dt 


+  Povs ' 


(Y-l)g  2 

— Vc0 


VP,s 


(19) 


Lo  ^ O 

Equations  15,  17,  and  19  are  a  complete  set  of  equations  for  the  propagation  of  infrasound  energy,  including  gravity 
and  wind,  in  a  form  amenable  to  FDTD  implementation. 


An  example  showing  the  sound  field  for  realistic  sound,  wind,  and  density  profiles  (Figure  2a)  is  shown  in  Figure 
2b,  for  a  source  at  100  km  altitude,  with  a  center  frequency  of  0.1  Hz.  Densities  and  temperature  profile  data  for  a 
location  near  the  I57US  infrasound  station  were  derived  from  the  MSIS-90  atmospheric  model,  made  available 
through  http://nssdc.gsfc.nasa.gov/space/  model/models/msis.html.  Wind  data  were  derived  from  the  Horizontal 
Wind  Model  (HWM)  1993  program  (Hedin  et.al.,  1996).  The  results  shown  in  Figure  3  are  for  the  pressure  field  at 
108  sec  after  the  "explosion."  As  indicated,  the  pressure  field  is  asymmetric  due  to  the  effects  of  wind,  although  the 
Mach  numbers  for  this  example  are  less  than  M=0.15  at  all  altitudes.  The  pressure  field  is  shown  in  the  panel  at  left, 
the  pressure  field  normalized  by  p0  "1/2  is  shown  in  the  panel  at  right.  As  indicated,  the  scaled  pressure  field  is 
reasonably  uniform  with  this  scaling. 


a) 


b) 


distanc  e  (km )  distant  a  ( km) 


Figure  2.  a)  Profiles  for  atmospheric  density,  sound  speed,  and  wind  speed,  b)  Pressure  field  (left)  and  the 

pressure  field  scaled  by  p0  1/2  (right)  for  a  source  with  center  frequency  at  0.1  Hz  and  at  an  altitude 
of  100  km,  for  medium  properties,  as  shown  in  Figure  2a. 


Linearized  acoustic  propagation,  including  attenuation,  wind  and  gravitational  effects 


It  has  been  shown  that  atmospheric  attenuation  is  approximately  proportional  to  the  square  of  frequency  (Sutherland 
and  Bass,  2004;  de  Groot-Hedlin,  2008a).  It  can  be  shown  (Batchelor,  1967),  that  the  effects  of  atmospheric 
attenuation  can  be  included  by  adding  two  terms  to  the  equation  for  conservation  of 


2_  1  _  2-  2_ 
momentum:  p(V  vs  +  —  V(V  •  vs)),  where  V  vsis  a  vector  with  components  V  v5 


:[V2VV2v 


s,yy2ysj’and 


^vs .  As  shown  in  de  Groot-Hedlin  (2008a),  the  first  term  yields  an  amplitude  decrease  proportional  to  the  square  of 
frequency  and  the  second  is  a  constant  as  a  function  of  frequency.  Thus,  the  equation  for  the  particle  velocity 
becomes 


dvo  _  _  _  __  _  o  1  _  _ 

Po(-Qf  +  vs  *  Vvo  +  vo  *  Vvs)  +  P5v0  •  Vv0  =  -Vps  - p sg  +  p(V  VS  +  -  V(V  •  vs))  +  ^vs  +  /  .  (20) 

Equations  20,  17,  and  19  are  a  complete  set  of  equations  for  propagation  of  infrasound  energy,  including  attenuation, 
gravity,  and  wind,  in  a  form  amenable  to  FDTD  implementation. 


An  example  comparing  infrasound  propagation  with  and  without  attenuation  through  a  model  with  constant  sound 
speed  and  exponentially  decreasing  density  (i.e.,  an  isothermal  atmosphere)  is  shown  in  Figure  3.  Attenuation  is 
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proportional  to  the  square  of  the  frequency  and  is  equivalent  to  that  at  an  altitude  of  1 10  km.  Pressure  waveforms  are 
sampled  at  locations  shown  above  and  below  the  source  and  are  compared  in  Figure  4,  along  with  associated 
waveforms. 


-0.05  0  0.05 


Range  (km) 


Range  (km) 


time  159  s 


o 

o 


o 


o 

5  10 


Range  (km) 


Range  (km) 


Range  (km)  Range  (km)  Range  (km) 


Range  (km)  Range  (km)  Range  (km) 


Figure  3.  (top)  No  attenuation.  Waveforms  are  sampled  at  receiver  locations  shown  by  the  circles.  The  source 
is  indicated  by  an  x.  (bottom)  same,  but  with  attenuation  equivalent  to  that  at  an  altitude  of  110  km. 


Figure  4.  (left)  Waveforms  are  shown  for  no  attenuation  (black),  constant  attenuation  (blue)  and  attenuation 
proportional  to  the  square  of  the  frequency  (red,  as  shown  in  Figure  3b).  (right)  associated  spectra, 
normalized  with  respect  to  peak  power  for  waveforms  propagated  through  non-attenuating  media. 
Colored  lines  become  lighter  with  increasing  distance  from  the  source. 
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Example  of  FDTD  applied  to  penetration  of  infrasound  into  shadow  zones 

Variations  in  wind  velocities  and  temperature  with  altitude  lead  to  the  refraction  of  acoustic  energy.  During  the  day, 
temperatures  and  adiabatic  sound  speed  typically  decrease  with  altitude,  causing  sound  to  refract  upward  (Reynolds, 
1873-1874).  The  effect  of  wind  on  refraction  is  often  equal  to  or  greater  than  that  of  temperature,  so  that  sound  is 
refracted  downward  in  the  direction  of  the  wind  and  upward,  away  from  the  ground,  in  the  opposite  direction. 


A  first-order  representation  of  acoustic  “zones  of  silence”  (ZoS)  can  be  obtained  using  the  high-frequency,  ray- 
theory  approximation  as  shown  in  Figure  5.  Rays  propagating  through  a  moving  atmosphere  are  computed  using 
equations  found  in  Garces  (1998).  Temperature  decreases  with  altitude  at  a  rate  of  6.5°C/km  (the  environmental 
lapse  rate)  and  the  wind  profile  was  chosen  to  approximately  offset  the  associated  decrease  in  sound  speed  with 
altitude.  Rays  were  launched  at  an  altitude  of  2  km  at  a  variable  interval  angle  such  that  the  ray  endpoints  would  fall 
at  0.5  km  increments  on  the  surface  for  a  homogeneous  velocity  model  (to  correct  for  cylindrical  spreading  at  the 
surface).  Therefore,  the  distance  between  ray  endpoints  in  Figure  5  (or  the  ray  density  per  unit  length)  is  a  measure 
of  the  pressure  amplitude,  and  shows  the  degree  of  focusing  and  defocusing  due  to  velocity  heterogeneity.  Pressure 
amplitude  is  slightly  greater  to  the  far  right  (in  the  wind  direction).  To  the  left,  pressure  amplitude  quickly  decreases 
starting  at  x  =  -5  km  and  is  effectively  undetectable  in  the  ZoS  beginning  at  about  x  =  -15  km.  In  nature,  the  first 
ZoS  away  from  a  near-surface  source  typically  extends  from  several  kilometers  from  the  source  to  between  120  and 
200  km,  depending  on  stratospheric  winds. 
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Figure  5.  A  “zone  of  silence”  (ZoS)  illustrated  by  the  high-frequency  approximation  of  ray  theory.  The 

assumed  adiabatic  sound  speed  (m/s)  and  horizontal  wind  speed  to  the  right  (m/s)  are  shown  in  (a) 
and  (b).  The  dashed  line  in  (b)  is  out-of-plane  wind  speed  (zero).  Every  tenth  ray  is  red.  Vertical 
exaggeration  is  6. 


Ray  theory  relies  on  a  high-frequency  approximation  to  wave  propagation.  So  its  use  assumes  that  sound  speeds 
vary  on  a  scale  length  much  larger  than  the  signal  wavelength.  This  approximation  starts  to  break  down  at 
infrasound  frequencies  under  typical  atmospheric  conditions,  where  sound  speed  gradients  of  4  m/s  per  km  of 
altitude  are  common.  A  more  complete  description  of  infrasound  propagation  near  the  ground  includes  finite- 
wavelength  effects  like  diffraction  (Embleton,  1996),  surface  waves  (Attenborough,  2002),  scattering  along  rough 
topography,  and  also  on  temporal  effects  like  turbulence  and  propagating  gravity  waves.  These  effects  lead  to  a 
penetration  of  acoustic  energy  into  the  ZoS  at  amplitudes  much  greater  than  that  predicted  by  ray  theory  (Embleton, 
1996;  Attenborough,  2002). 

Zones  of  silence  (ZoS)  can  be  examined  in  better  detail  using  numerical  modeling  algorithms  that  take  dispersion, 
scattering,  and  other  linear,  full-wave  effects  into  consideration.  For  instance,  snapshots  of  the  evolving  pressure 
field  are  shown  in  Figure  6  from  FD  modeling  using  the  same  velocity  model  shown  in  Figure  5.  The  source  (again 
at  a  2-km  altitude)  consists  of  the  superposition  of  a  0.9  Hz  and  a  0.2  Hz  wavelet,  multiplied  by  a  cosine  taper. 
Waveforms  are  shown  at  5-km  increments  along  the  ground  from  -30  km  to  +30  km  (Figure  7).  The  waveforms  are 
scaled  by  the  square  root  of  the  source-receiver  distance  to  correct  for  spreading.  Peak  amplitudes  are  generally 
equal  to  the  right,  in  agreement  with  the  prediction  from  ray  theory.  However,  the  decrease  in  peak  amplitude  to  the 
left  is  not  as  severe  within  the  ZoS  as  predicted  by  ray  theory.  Lower  frequencies  are  diffracted  into  the  ZoS;  the 
power  spectra  for  each  waveform  show  that  0.9  Hz  energy  decreases  more  rapidly  than  that  around  0.2  Hz  within  the 
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ZoS.  Therefore,  ZoS  predicted  by  ray  theory  may  be  completely  penetrated  by  relatively  lower- frequency 
infrasound,  or  exhibit  frequency  dependence  with  penetration  distance.  This  is  in  line  with  observations:  numerous 
studies  have  observed  infrasound  energy  in  classical  shadow  zones  (e.g.,  Hedlin  et  al.,  2008b;  Ottemoller  and  Evers, 
2008). 
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Figure  6.  Pressure  snapshots  at  38.4,  57.6  and  76.8  seconds,  derived  from  FD  numerical  modeling.  The  wind 
and  sound  speeds  are  the  same  as  for  Figure  4.  The  source  is  at  2  km  in  altitude,  marked  by  a  red 
asterisk,  and  receiver  locations  are  marked  by  an  x. 


Figure  7.  (top)  Waveforms  predicted  from  FD  modeling  for  a  source  near  a  flat  surface.  Waveforms  are 
corrected  for  spreading  at  receivers  located  at  5-km  intervals  along  the  ground,  (bottom)  Power 
spectra  are  normalized  by  the  level  of  0.2  Hz  energy. 
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The  FD  algorithm  was  re-run,  this  time  with  variable  topography  to  determine  its  effect  on  infrasound  propagation 
into  a  shadow  zone.  The  ground  surface  was  modeled  as  a  sine  wave  with  a  wavelength  of  10  km,  and  a  maximum 
height  of  500  m.  Receivers  are  located  along  the  air/ground  interface  at  an  elevation  of  z=0  m,  since  the  receivers  all 
lie  along  the  zero  crossings  of  the  sine  wave.  Several  snapshots  are  shown  in  Figure  8  and  the  waveforms  and 
spectra  are  shown  in  Figure  9.  As  shown,  the  introduction  of  topography  complicates  the  response.  The  waveforms 
are  amplified  at  some  point  (e.g.,  at  -10  km  and  at  +15  km)  and  diminished  at  others  (-15  km  and  +10  km  and  +20 
km).  The  higher  frequencies  are  diminished  more  than  the  low  frequencies.  The  results  suggest  that  ZoS  depends  not 
only  on  atmospheric  temperatures  (hence,  adiabatic  sound  speed)  and  temperatures,  but  is  also  a  complicated 
function  of  frequency  and  topography.  The  addition  of  topography  appears  to  further  strip  off  high  frequencies  in 
comparison  to  the  case  for  flat  topography  (Figures  6  and  7  above). 
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Figure  8.  Pressure  snapshots  at  38.4  and  57.6  seconds,  derived  from  FD  numerical  modeling  for  the 

air/ground  interface  modeled  as  a  sine  wave.  The  wind  and  sound  speeds  are  the  same  as  for  Figure 
4.  The  source  is  at  2  km  in  altitude,  marked  by  a  red  asterisk;  receiver  locations  are  marked  by  x. 
Receivers  are  at  0-km  altitude,  as  in  the  previous  example. 


Figure  9.  (top)  Waveforms  predicted  from  FD  modeling  for  a  source  near  a  sinusoidally  varying  surface. 

Waveforms  are  corrected  for  spreading  at  receivers  spaced  at  5-km  intervals  along  the  ground, 
(bottom)  Power  spectra  are  normalized  by  the  level  of  0.2  Hz  energy. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


An  algorithm  has  been  developed  to  allow  for  the  combined  effects  of  wind  and  gravity  on  the  infrasound  field.  The 
addition  of  attenuation  into  this  code  will  be  accomplished  soon.  The  program  is  currently  running  slowly  on  a 
matlab  platform.  The  plan  is  to  convert  it  to  fortran  in  the  near  future  for  considerable  improvements  in 
computational  speed.  After  that,  the  code  will  be  extended  to  three  dimensions  and  tested  against  real  data. 
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